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Search for the decay B + —> K*°K* + at Belle 
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We report a search for the rare charmless decay B + —> K*°K* + using a data sample of 772 x 10 6 
BB pairs collected at the T(45) resonance with the Belle detector at the KEKB asymmetric-energy 
e + e“ collider. No statistically significant signal is found and a 90% confidence-level upper limit is 
set on the decay branching fraction as B(B + —> K*°K * + ) < 1.31 x 10 -6 . 


seem to follow the expected pattern [4] . The polarization 
puzzle is a prime motivation for measurements in other 
B —► VV decays to test predictions of the QCD factor¬ 
ization and perturbative QCD approach. The sensitivity 
to f L is obtained by considering the decay process in the 
helicity basis. In the B + —i K*°K* + decay, this basis is 
defined with the two K* rest frames in which the helicity 
angles d K ,+ and d K *o are measured between the daugh¬ 
ter momentum (K ± or ti ) of each I\* and the direction 
opposite the B meson. 

The B + —> K*°K* + decay proceeds via electroweak 
and gluonic b —>■ d, loops, as shown Fig. [Tj The expected 
branching fractions for B meson decays to VV final states 
are calculated in several papers [SHU]. The branching 
fraction of B + —> K*°K* + is predicted to be (0.1 — 1.1) x 
10~ 6 in QCD factorization [BJ ECO and (0.3 — 0.9) x 10~ 6 
in perturbative QCD 

The BABAR Collaboration has measured the longitu¬ 
dinal fraction f L = 0.75t(y26 ±0.03 and the branching 
fraction B = (1.2 ± 0.5 ± 0.1) x KT 6 for B+ -> K*°K*+ 
using a data sample of 467 x 10 6 BB pairs [T2j, where 
the first uncertainty is statistical and the second is sys- 


PACS numbers: 13.25.Hw, 11.30.Er, 12.15.Hh 

The study of charmless B meson decays provides 
a powerful probe to search for new physics [T] be¬ 
yond the standard model. We search for B + —> 
RT*°(892)AT* + (892), a B —> VV decay channel medi¬ 
ated by the b — i d transition for which the so-called 
polarization puzzle is yet to be solved; here, V denotes 
a vector meson. A naive counting rule for light vector 
mesons predicts the longitudinal-polarization fraction to 
be /i ~ 1 — Oirriy/m 2 B ) in such decays [2j. However, 
in loop-dominated modes such as B —> cf>K* 0, the f L 
values are found to differ significantly from this predic¬ 
tion. In contrast, tree-dominated decays, e.g., B —► pp 
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FIG. 1: (a) Electroweak and (b) gluonic b —¥ d penguin loop 

diagrams for B+ —> 
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tematic. It has also obtained the A 0 —> X*°K*° decay 
branching fraction B = ( 1 . 28 ^q' 3 q ± 0.11) x 10~ 6 113] . 
On the other hand, Belle reported an upper limit on the 
branching fraction for B° -> K*°K*° (A 0 -> K*°K*°) of 
0.81 x 10 -6 (0.20 x 10 -6 ) ill]. Owing to the smallness 
of the underlying CKM matrix elements, the b —► d tran¬ 
sitions (dominant in B —> K*K* decays) are suppressed 
compared to b —> s and hence the related channels are 
not so well measured. Therefore, precise measurements 
based on high statistics are needed to shed more light on 
the polarization puzzle. 

Our results are based on a data sample containing 
772 x 10 6 BB pairs, corresponding to an integrated lu¬ 
minosity of 711 fb —1 , recorded at the T(45) resonance 
with the Belle detector m at the KEKB asymmetric 
energy e + e _ (3.5 on 8.0 GeV) collider ITS]. The princi¬ 
pal detector components used in the study are a silicon 
vertex detector (SVD), a 50-layer central drift chamber 
(CDC), an array of aerogel threshold Cherenkov coun¬ 
ters (ACC), a barrel-like arrangement of time-of-flight 
scintillation counters (TOF), and a CsI(Tl) crystal elec¬ 
tromagnetic calorimeter (ECL). All these components are 
located inside a 1.5 T solenoidal magnetic field. Two in¬ 
ner detector configurations were used: a 2.0 cm beampipe 
and a 3-layer SVD for the first sample of 152 x 10 6 BB 
pairs, while a 1.5 cm beampipe, a 4-layer SVD and a 
small-cell CDC to record the remaining 620 x 10 6 BB 
events m- The latter sample has been reprocessed with 
an improved track reconstruction algorithm, which sig¬ 
nificantly increased the signal detection efficiency. 

The B + —> k*°K* + candidate is reconstructed from 
the subsequent decay channels of AT* 0 —> K~tt + and 
K* + —> K + tt° (A'g7r + ), where K* refers to the AT* (892) 
meson [IHl- 

Charged tracks are required to have a transverse mo¬ 
mentum greater than 0.1 GeV/c and an impact parame¬ 
ter with respect to the interaction point less than 0.3 cm 
in the j— (f) plane and 4.0 cm along the 2 axis. Here, the 
2 axis is the direction opposite the e + beam. Charged 
kaons and pions are identified by means of a likelihood 
ratio Rk/-k = £k/(£k + where £/<(£ w ) denotes 
the likelihood for a track being due to a kaon (pion). 
These likelihoods are calculated using specific ionization 
in the CDC, information from the TOF, and the number 
of photoelectrons from the ACC. Kaon identification ef¬ 
ficiencies are 98.1% (99.0%) for transversely and 97.2% 
(97.5%) for longitudinally polarized cases, and pion iden¬ 
tification efficiencies are 97.2% (98.6%) for transversely 
and 97.3% (98.9%) for longitudinally polarized cases in 
the AT* + — > ATg7r + (AT* + — > K + ir °) channel. Fake rates 
for kaons and pions are approximately 0.1% and 0.8%, 
respectively. 

Neutral n° and Kg mesons are reconstructed with a 
pair of photons and charged pions, respectively. The 7 r° 
candidates are required to have each daughter photon’s 
energy greater than 0.05 GeV (0.10 GeV) for the barrel 


(endcap) region of the ECL, a reconstructed invariant 
mass in the range 0.118 GeV/c 2 < m 77 < 0.150 GeV/c 2 , 
and a 7r° mass-constrained fit statistic, y; 2 0 , smaller than 
50. The mass requirement corresponds to ±3er around 
the nominal 7r° mass |TS] ■ The Kg candidates are selected 
with the following criteria. The 2 -distance between the 
two helices at the 7r + 7r - vertex position must be less than 
2.5 cm. After this initial selection, the pion momenta 
are refitted with a common vertex constraint. The flight 
length of the Kg candidate must lie between 2 and 20 cm. 
The impact parameter with respect to the interaction 
point must be greater than 0.1 cm in the r — (f> plane. 
Finally, we require the reconstructed invariant mass to 
be in the range 0.478 GeV/c 2 < m nn < 0.516 GeV/c 2 , 
corresponding to ±5cr around the nominal K° mass m- 

The K* candidates are reconstructed by defining the 
mass range from 0.78 to 1.00 GeV/c 2 that corresponds to 
approximately ±2.1cr around the nominal K* mass Il9l . 
In order to reduce the contribution of misreconstructed 
candidates in the AT* + —► K + n n decay, we require the 
helicity angle of the AT* + candidate to satisfy cos e K ,+ < 
0 . 8 . 

We define two kinematic observables in the form of 
the energy difference (A E = As — Abeam) and the beam- 
energy constrained mass (M bc = JrV^beam “ IAb| 2 c 2 ), 
where A be am and Eg ( Pb ) are the beam energy and the 
energy (momentum) of the B meson candidate, respec¬ 
tively, in the e + e~ center-of-mass (CM) frame. For the 
K*+ — ). K + n° channel, where the AA resolution is poor 
due to shower leakage in the ECL (20], we use the follow¬ 
ing quantity instead of M bc : 


( -* _\ 2 2 

Py.oC + |pT + | \J (ffbeam - %o) 2 - 

where mjf.+ is the K* + mass. We retain B candidates 
that satisfy |AA| < 0.15 GeV and M^} > 5.25 GeV/c 2 . 

The dominant background arises from the e + e _ —> 
qq (q = u, d, s, c) continuum process. To suppress these 
events, a neural network m is employed by combin¬ 
ing the following four quantities: a Fisher discriminant 
formed from 16 modified Fox-Wolfram moments 122] . the 
cosine of the angle between the momentum of signal B 
candidate and the 2 axis in the CM frame, the separa¬ 
tion along the 2 axis between the vertex of the signal B 
and that of the recoil A, and the recoil A’s flavor-tagging 
information [25j. To reconstruct the decay vertex of the 
recoil A, the tracks not associated with the signal A are 
used. The training and optimization of the neural net¬ 
work are accomplished with signal and continuum Monte 
Carlo (MC) simulated events. The signal MC sample is 
generated with the EvtGen program [M], taking final- 
state radiation effects into account via PHOTOS . 125] . 
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The neural network output ( Cnb ) ranges from —1 to 
+1; an event near +1 (—1) is more signal (continuum)- 
like. We require Cnb > —0.5 to reduce substantially the 
amount of continuum background. This requirement pre¬ 
serves approximately 94.7% (94.5%) of the signal while 
suppressing 75.6% (71.2%) of the continuum background 
in K* + —► K + tt° (K*+ —> Kgn + ). As the remainder 
of the Cnb distribution has a sharp peak near unity, we 
use a transformed quantity to enable its modeling with 
an analytic shape: 


C'nb ~ 1 °§ 


NB ^NB \ 

C m ax I ’ 

NB ^NBJ 


( 2 ) 


where C$§ = -0.5 and C%% = 0.997 (0.995) in K*+ -> 
K+tt° (K*+ K%ir+). 

After all selection criteria are applied to the signal 
MC sample, the average number of signal candidates 
per event is 1.16 (1.13) for longitudinally (transversely) 
polarized decays in I\* + —> K + ir° and 1.10 (1.06) in 
K* + —» KgTT + . We choose the candidate having the 
smallest x 2 0 + X% ( X 2 K o + Xb) value in K * + K+7T ° 

( K* + —> K^n + ), where the B vertex is obtained by 
charged tracks except for those from Kg and (x'k° ) 

is the B (Kg) vertex-fit statistic. We refer to the right- 
combination (RC) as the correctly reconstructed B me¬ 
son decays and the self-crossfeed (SCF) as the misrecon- 
structed signal component. MC simulations show that 
the SCF fraction is 15.5% (10.2%) for the longitudi¬ 
nally (transversely) polarized case in K* + —> K + n° and 
7.7% (3.5%) for the longitudinally (transversely) polar¬ 
ized K* + —> KgTT + decay. 

The charm BB background originating from the b —> c 
transition remains after all event selection criteria are ap¬ 
plied. In the MC sample, we find no peaking structure in 
A E, MA , and the invariant masses formed by combining 
two or three final-state particles. We also do not observe 
any specific charm decay mode in this sample. The other 
possible backgrounds are largely due to b —> u, d, s tran¬ 
sitions from charmless B decays. This background has 
no peaking structure in the signal enhanced region of 
|AA| < 0.05 GeV, while a peaking structure originated 
from B + —> p°K* + and B + —> ttttK* + with K* + —> 
Kgn + is seen at A E ~ 0.07 GeV. Other backgrounds 
involving higher K* states such as K*K% (1430) and 
A*R”q(1430), KnK* decays, and the nonresonant four- 
body KnEgir ( KttKtt °) decays also contribute. The 
A'*A'|(1430) decays are simulated based on the theoret¬ 
ical expectations [26] for branching fractions and polar¬ 
izations. The contributions of A'*Aq( 1430) decays are 
estimated on both K* mass sidebands, where the Kgn 
(Kn°) mass sideband is 0.78 GeV/c 2 < ttik-k < 1-00 
GeV/c 2 and 1.00 GeV/c 2 < m K o^^o) < 1.52 GeV/c 2 
and the I\n mass sideband is 1.00 GeV/c 2 < mx-n < 1-52 
GeV/c 2 and 0.78 GeV/c 2 < m K ^(Kir 0 ) < 1-00 GeV/c 2 . 
The B + —>■ c/)K* + background arising from pion-to-kaon 


misidentihcation is suppressed by rejecting events with 
an invariant mass of the K + K~ pair between 1006.5 and 
1032.5 MeV/c 2 . 

The B —>• VV decay rate does not depend strongly on 
the azimuthal angle, (f >, between the two decay planes of 
the vector mesons. Therefore, it can be integrated out to 
obtain the differential decay rate [25] 


1 ci 2 r 
T dcos9-^, 0 dcos8 K *+ 


^j(l - f L )sin 2 8 1? ,oSm 2 e K ,+ 

9 

+ -f L cos 2 9j ? , 0 cos 2 9 K ,+ . (3) 


We obtain the branching fraction B and the longitu¬ 
dinal polarization fraction using a simultaneous fit to 
the K* + —> Kgn + and K* + —> K + n° decay channels. 
This is an unbinned extended maximum likelihood (ML) 
fit to the distributions of A E and M h J , the invariant 
mass and the cosine of the helicity angle of the two K* 
candidates, and C' NB . The extended ML function for 
each decay channel is 




N 

n 

i =1 




a, 


(4) 


where Vj(xi\dj) is the product of uncorrelated one- 
dimensional (ID) probability density functions (PDFs) 
for event category j, calculated for the seven mea¬ 
sured observables Xi of the i-th event, rij is the yield 
for this event category, and N is the total number of 
events. The parameters dj describe the expected dis¬ 
tributions of the measured observables for event cate¬ 
gory j, and are extracted from MC simulations and the 
(A* mass, M^*J) sideband data. For the simultaneous 
fit, the total likelihood is obtained by multiplying the 
likelihoods for the K* + —> Kgir + and K* + —> A' + 7r° 
decay channels (indexed by k). With an assumption 
of equal production of B + B~ and B°B° pairs at the 
T(4S') resonance, the signal yield of channel k is given 


by n sig , k = BX f L ^rec,k + i 1 ~ f l) 


c rec ,k 


x n B k x n bb , 


where N BB is the number of BB pairs, n s i g is the num¬ 
ber of signal events, and IL B k is the product of the sub¬ 
branching fractions. The detection efficiency for the Ion- 
gitudinally (transversely) polarized mode, e re c is equal 
to 11.58 ± 0.02% (14.41 ± 0.02%) and 12.35 ± 0.02% 
(17.29±0.02%) for the K*+ -> K%n+ and K*+ -> K+n° 
channels, respectively. These are determined primar¬ 
ily from the signal MC sample and then corrected for 
a modest difference of kaon-identification efficiency be¬ 
tween data and simulations, given by r K /n = s < x/^/s^9 n , 
where (e^^.) is the efficiency of the Rx/-k require¬ 
ment in data (simulations). The r^/, value per charged 
pion (kaon) track is 0.96 (1.00), resulting in a total effi¬ 
ciency of 0.92 (0.96) for K* + —>■ Kgir + (K* + —> K + ir°). 
Though mild linear correlations of up to 15% exist in 
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TABLE I: List of PDFs used to model the AE, M^* c \ m K° 7r(K+ir 0 )^ cos ^Km cosOK s^{k+-k°) an d ^nb distributions for the 

various event categories. G, AG, CB, ARG, (r)BW, P i, LASS, Hist and Erf stand for Gaussian, asymmetric Gaussian, Crystal Ball m, 
ARGUS function |28| . (relativistic) Breit-Wigner function, z-th order Chebyshev polynomial, LASS parameterization for the Aq( 1430) 
line shape, histogram and error function, respectively. Two different PDFs are used to model cosQk-k on the two samples of tyik-k < 0.83 
GeV/c 2 and m Kl r > 0.83 GeV/c 2 . 


Final state 

Event category 

AE 


m K TT 

m K°Tr(K + nO) 

COS0KTT 

COS @K°Tr(K+TT°) 

C'nb 


Signal (RC) 

2G 

G 

rBW 

rBW 

Hist / Hist 

Hist 

2AG 


Signal (SCF) 

Hist 

Hist 

Hist 

Hist 

Hist / Hist 

Hist 

AG 


Continuum qq 

Pi 

ARG 

rBW+Pi 

rBW+Pi 

P 6 x Erf / P 4 x Erf 

PgxErf 

2G 


Charm BB 

Pi 

ARG 

Pi 

P 2 

P 4 / P 4 xErf 

Ps 

AG 

AT - 7T + Ag7r + 

Charmless BB 

C+P 2 

G+ARG 

BW+Pi 

BW+Pi 

Hist / Hist 

Hist 

AG 


B+ -+ (K^)o°K*+ 

2G 

2G 

LASS 

rBW 

Hist / Hist 

Hist 

AG 


B+ -+ KfK* + 

2G 

2G 

BW 

rBW 

Hist / Hist 

Hist 

AG 


B + -+K*°(K-k)*+ 

g+p 2 

G 

rBW 

LASS 

Hist / Hist 

Hist 

AG 


B + K*°K* + 

g+p 2 

G 

rBW 

BW 

Hist / Hist 

Hist 

AG 


B + —y four-body 

g+p 2 

g+p 2 

Pi 

Pi 

Hist / Hist 

Hist 

AG 


Signal (RC) 

CB+G 

CB 

rBW 

rBW 

Hist / Hist 

Hist 

2AG 


Signal (SCF) 

Hist 

Hist 

Hist 

Hist 

Hist / Hist 

Hist 

AG 


Continuum qq 

Pi 

ARG 

rBW+Pi 

rBW+Pi 

P 6 xErf / P s xErf 

Pe 

2AG 


Charm BB 

P 2 

ARG 

Pi 

Pi 

P 4 / P 4 xErf 

P 4 

AG 

K-n+K+n 0 

Charmless BB 

P 4 

P 4 

BW+Pj 

BW+Pi 

Hist / Hist 

Hist 

AG 


B+ -+ (K^)* 0 °K*+ 

cb+p 2 

2G 

LASS 

rBW 

Hist / Hist 

Hist 

2AG 


B+ -+ KfK*+ 

cb+p 2 

2G 

BW 

rBW 

Hist / Hist 

Hist 

2AG 


B + —>7r ,0 (RV)*+ 

cb+p 2 

2G 

rBW 

LASS 

Hist / Hist 

Hist 

2AG 


B + k*°k;+ 

cb+p 2 

2G 

rBW 

BW 

Hist / Hist 

Hist 

2AG 


B + —¥ four-body 

g+p 2 

g+p 2 

Pi 

Pi 

Hist / Hist 

Hist 

AG 


the signal, such as between (A E, Mb c ), their contribu¬ 
tions to the fit bias (described later) due to our use of 
uncorrelated ID PDFs are negligible. 

Table [i] lists the PDF shapes used to model AE, M^*J, 
mKiT , %»!(&«)) cos 6 K tt, cos d Ks ^ K ^ and C' NB for 
different event categories. We fix the parameters of the 
RC signal PDF shapes to the MC values. We correct the 
parameters of the RC signal AE, M^*j and C' NB PDFs to 
account for modest data-MC differences; the correction 
factors are obtained from a high-statistics control sample 
of B + —» [i~)K* + . The same calibration factors 

are also applied to the higher- K* and nonresonant back¬ 
grounds. 

The continuum background PDF parameters that are 
allowed to vary are the slope of AE , the shape of ME, 
the fraction of the relativistic Breit-Wigner function, the 
polynomial coefficients of the K* masses, and the mean 
and two widths for the core asymmetric Gaussian func¬ 
tion of C' NB . All other PDF parameters are fixed and 
determined from MC samples. We use an error function 
to describe the falling reconstruction efficiency due to 
low-momentum tracks in the continuum as well as the 
BB helicity angle distributions. We use the simulta¬ 
neous fit with two different cos 9ktt PDFs, correspond¬ 
ing to the two samples of tuk-k < 0.83 GeV/c 2 and 
iriKn > 0.83 GeV/c 2 , to treat the correlation between 
rriKir and cosOk-k that originates from the B —► <j)K* 


veto. 

The yields for all event categories except for the rel¬ 
ative amount of SCF to RC signal, the charmless BB, 
higher K* and nonresonant background components are 
allowed to vary in the fit. We fix the yields of charmless 
BB backgrounds based on a high-statistics MC sample, 
which includes possible charmless rare B decays. In or¬ 
der to validate our fitting procedure, we perform the fit to 
ensembles of 500 pseudoexperiments using the extracted 
fitted yields from data and events of all components that 
are arbitrarily chosen from the simulated MC samples. 
We obtain and correct for fit biases of 1.8% and 8.2% for 
B and /l, respectively, and assign 50% of each bias as its 
systematic uncertainty. 

The A)) (1430) resonance, together with an effective- 
range nonresonant component, are modeled with the 
LASS function [30 !. whose parameters are taken from 
Ref. HQ. Yields of (. K^)* 0 °K*+, K*°{Kn)* 0 + and 
four-body decay backgrounds are measured by a simul¬ 
taneous fit to the sidebands of the two K* masses. 
To combine the results of the two K* + decay chan¬ 
nels, both fits share the branching fraction parameters 
of (A'7r)g°A'* + , A'*°(A'7r)g + and four-body decay back¬ 
grounds for K* + — > K b tt + and AT* + — > K + tt° in the 
simultaneous fit. In the fit, these background yields in 
the K* mass signal region from 0.78 to 1.00 GeV/c 2 are 
estimated from the K* mass PDFs on the two K* mass 
sidebands. 
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FIG. 2: (color online). Projections for B + —> K *°(— >■ K n+)K* + (—>■ i^^ 7 r + ) of the multidimensional fit onto A E, Mb c , K*° mass, 
K*+ mass, cosine of K*° helicity angle, cosine of K*+ helicity angle, and C' NB for events selected in a signal enhanced region with the 
plotted variable excluded. Points with error bars are the data, the solid curves represent the full fit function, the hatched regions are the 
signal, the dashed curves show the combined continuum and BB backgrounds, and the dotted curves are the higher K* and nonresonant 
backgrounds. 


We obtain the biases of the (Ktt)qK* and four-body 
decay yields by applying the fit to ensembles of 500 
pseudoexperiments using the extracted fitted yields from 
the K* mass sidebands. Fit biases for the yields of 
K * + , K*°(Kn)* 0 + and four-body decays are, re¬ 
spectively, 3.0 (2.6), 2.0 (2.0) and 0.8 (0.4) in the K* + —y 
Kgir + ( K* + —> K+tt°) sample. We correct for the fit 
biases and assign 50% of each to the systematic un¬ 
certainties. The measured yields in the K* mass side¬ 
bands are extrapolated to the K* mass signal region us¬ 
ing the K* mass PDFs. We obtain the background yields 
N(k^)*o k * + = 1 - 9 ^ 2.8 ^K*°{Kir)l+ = 3 - 3 - 2.3 

(3.2 ± 1.9), and iV^body = 2.5 ± 3.0 (1.2 ± 1.4) in the 
K*+ KgTr + ( K* + — > K + tt°) samples, where errors 
are a quadratic sum of the statistical and systematic un¬ 
certainties. 

The total event sample for B + —> K*°K* + consists 
of 23338 and 50212 events with K* + —> K B n + and 
K*+ —> K + n°, respectively. The result of the ML fit is 
summarized in Table [TTJ We take the sub-branching frac¬ 
tions B(K*° -± K~ir+) = 2/3, B(K*+ -+ K 0 n+) = 2/3, 
B{K*+ -± K+ir°) = 1/3 and B{K° -> K% -> tt+tt") = 
0.5 x (69.20 ± 0.05)% [19]. The signal significance S is 
dehned as v /-21og(£ max /£ 0 ), where £ max (£ 0 ) is the 
likelihood value when the signal yield is set to its nomi¬ 
nal value (zero). The systematic uncertainty (discussed 
below) is included in this significance calculation by con¬ 
volving the statistical likelihood with an asymmetric 
Gaussian distribution whose width equals the total sys¬ 
tematic error. The total significance of the signal yield is 
2.7 standard deviations (a). The upper limit (UL) on the 
branching fraction is calculated at 90% confidence-level 
(CL) by using the formula / 0 8ul C(B)dB/ J 0 °° C(B)dB = 


0.9. The result is Sul = 1.31 x 10 6 . 

Figures [2] and [3] show the projections of the two fits 

— (*) —' 

onto AE , ME, K* masses, cosines of the helicity angle 
and C' NB for I\* + — ► Ks n+ aR d K * + K + n°. The 

candidates and PDFs in each figure are projected in the 


TABLE II: Summary of results for the fitted yields, average effi¬ 
ciencies e rec for the fitted /L, sub-branching fractions lon¬ 

gitudinal polarization fraction fi,, branching fraction J3(B~^ —»■ 
K*°K*+), signal significance S , and B upper limit at 90% CL. 
The first error is statistical and the systematic error is quoted last, 
if given. 


Final state 

+ 

coc 

+ 

K-tt+K+tt 0 

Yields (events): 

Total 

23338 

50212 

Signal 

15.81 li 

16.7/™ 

qq 

22982/^ 

49733/2to 

Charm BB 

265/(49 

290/(g® 

Charmless BB (fixed) 

78 

166 

(K^y 0 °K* + (fixed) 

1.9 

1.6 

K*°{Ktt)* 0 + (fixed) 

3.3 

3.2 

KfK*+ (fixed) 

0.45 

0.30 

K*°KZ+ (fixed) 

0.10 

0.06 

four-body decay 

2.5 

1.2 

Efficiencies: 

£rec(%) 

11.58 ±0.02 

12.35 ±0.02 

n&(%) 

15.37 

21.96 

Results: 



h 

1.06 ±0.30 ±0.14 

B(xl0" 6 ) 

0.77l°;3o ± °- 12 

S(a) 

2.7 

13(xl0~ 6 ) upper limit (90% C.L.) 

1.31 
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FIG. 3: (color online). Projections for B + — > K*°( -+ K 7r + )K*+(— * K + tt°) of the multidimensional fit onto A E, M£ c , K*° mass, K* + 
mass, cosine of K*° helicity angle, cosine of K* + helicity angle, and C' NB . The same projection criteria and legend as Fig. [ 2 ] are used. 


signal-enhanced region: |A.E| < 0.05 GeV, > 5.27 
GeV/c 2 , 0.83 GeV/c 2 < mx* < 0.95 GeV/c 2 &ndC' NB > 
3. 

One of the sources for fit bias is its inaccurate esti¬ 
mation (based on the ensemble test) due to the limited 
size of the qq MC samples. The MC samples are gener¬ 
ated under a partially correlated qq PDF. We estimate 
an additional fit bias from the results of a comparison 
between ensemble tests using limited and correlated MC 
samples. We calculate the total fit bias uncertainty as 
the quadratic sum with this additional fit bias. The un¬ 
certainties due to the fixed yields for the higher K* and 
nonresonant backgrounds are estimated by varying the 
corresponding yields by their errors. The uncertainties 


TABLE III: Summary of systematic uncertainties (%) on the 
branching fraction and longitudinal polarization fraction. 



B 

fL 

Fit bias 

4.72 

6.81 

PDF modeling 

5.40 

5.32 

Histogram PDF 

1.07 

1.11 

Calibration factors 

+5.32 

-3.82 

- 

Track reconstruction 

2.10 

- 

PID efficiency 

3.90 

- 

Cnb efficiency 

1.08 

- 

Kg reconstruction 

0.73 

- 

7T° reconstruction 

4.09 

- 

Fractions of misreconstructed events 

+3.32 

-1.48 

+ 1.92 
-2.80 

Nonresonant & higher K* background 

+9.54 

-9.73 

+3.76 

-4.10 

Limited MC statistics 

0.31 

- 

Charmless BB background 

+2.13 

-0.67 

- 

Number of BB events 

1.37 

- 

Interference with (Kn)o 

5.80 

9.69 

Total 

+ 16.2 

-15.4 

+ 13.7 
-13.9 


due to the fixed fractions of misreconstructed events are 
calculated by varying them by ±50%. The charmless BB 
background yield is also varied by a conservative ±50% 
to cover any mismodeling of these processes in the MC 
sample used to estimate the yield. The change in the 
signal yield is taken as the systematic uncertainty. 

We estimate the effect of possible interference be¬ 
tween the K* and spin-0 final states [nonresonant and 
A'g(1430)] by including interference terms with variable 
phases in the relativistic Breit-Wigner function of the 
spin-0 final-state mass. In this estimation, we assume 
the K* helicity angle distributions for /i = 0 and /l = 1 
in the KnK* decay to be the same as those of our signal 
decay. We vary the amplitude and phase of the inter¬ 
ference term and the fractions of + = 0 and f L = 1 
components of KnK* from 0 to 1. We assign the result¬ 
ing shifts as the systematic uncertainties after refitting 
with this modified function. 

We vary the bin height for all histogram PDFs by its 
statistical error and repeat the fit. Deviations from the 
nominal fit are added in quadrature to determine the 
uncertainty due to the histogram PDFs. The PDF mod¬ 
eling uncertainty is obtained by varying the fixed shape 
parameters by their errors. We assign an uncertainty on 
the absolute scale of the reconstruction efficiency due to 
the limited signal MC statistics. The uncertainty due to 
calibration factors to correct for the difference between 
data and simulations is obtained by varying those factors 
by their errors. We assign an uncertainty due to the dif¬ 
ferent continuum suppression efficiencies at Cnb = —0.5 
in data and MC by using the B + —> J/+++ ^jl~)K* + 
control sample. We also include reconstruction efficiency 
uncertainties for charged tracks (0.35% per track) by 
using partially reconstructed D* + —> D°(Kgn + n~)n + , 
particle identification (PID) uncertainties by using the 














































































D*+ —» D°(K~tt + )tt + control sample, and the uncer¬ 
tainty on the number of BB pairs. The systematic un¬ 
certainty due to the n° reconstruction is obtained by 
comparing data-MC differences of the yield ratio between 
—> TT°n°Tr 0 and ij — > n + Tr~n 0 . The systematic uncer¬ 
tainties on the branching fraction and longitudinal polar¬ 
ization are listed in Table Hill 

In summary, we have searched for the charmless 
hadronic decay B + —► K*°I\* + using the full BB pair 
sample collected with Belle. We find a 2.7cr excess 
of signal with a branching fraction B = (0.77+ 0 ’ 3 q ± 
0 .12) x 10 1 2 3 4 5 6 7 and a longitudinal polarization fraction 
fh = 1.06 ± 0.30 ± 0.14. We obtain a branching frac¬ 
tion upper limit of 1.31 x 10~ 6 at 90% CL. 
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